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Abstract 

We calculate analytically asymptotic values of quasi-normal frequencies of four- dimensional 
Kerr black holes by solving the Teukolsky wave equation. We obtain an expression for arbitrary 
spin of the wave in agreement with Hod's proposal which is based on Bohr's correspondence 
principle. However, the range of frequencies is bounded from above by 1/a, where a is the 
angular momentum per unit mass of the black hole. Our argument is only valid in the small-a 
limit which includes the Schwarzschild case. 
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Quasi-normal modes of black holes in asymptotically flat space-times play an important role 
in black hole physics and have attracted a lot of attention recently [1-11]- For high overtones, the 
imaginary part of the quasi-normal frequencies may be easily understood in terms of the poles 
of thermal Green functions, the spacing of frequencies being 27riTH, where Th is the Hawking 
temperature. On the other hand, the real part approaches an interesting asymptotic value 

Uu^ThIuS (1) 

in the case of gravitational perturbations. This has long been known numerically [12-16]. The 
analytical expression (1) was proposed by Hod [17] and was subsequently shown to be related 
to the Barbero-Immirzi parameter [18, 19] of Loop Quantum Gravity (see [20-24] and references 
therein). The starting point of these arguments is the entropy-area relation S = ^ A [25], 
relating the number of microstates of the black hole, S, to the area of the horizon, A. Mukhanov 
and Bekenstein [26] proposed that the area spectrum of black holes be discrete with spacing of 
eigenvalues 

5A^AG\nk , A; = 2,3,... (2) 

in units such that h = c = 1. Hod [17] used Bohr's correspondence principle to relate the real 
part of quasi-normal frequencies (1) to the area spectrum (2). His argument suggested that 
k — 3 in (2) instead of the expected k — 2 [26] . This is intriguing from the loop quantum gravity 
point of view because it suggests that the gauge group should be SO (3) rather than SU{2). It 
seems that quasi-normal modes may lead to a deeper understanding of black holes and quantum 
gravity. 

The asymptotic expression (1) has been derived analytically by Motl and Neitzke [3, 5] who 
used an interesting monodromy argument which relied heavily on the unobservable black hole 
singularity. First-order corrections have also been calculated analytically through a WKB analy- 
sis [8] for a gravitational wave and by solving the wave equation perturbatively [27] for arbitrary 
spin of the wave. 

Extending the above results to rotating (Kerr) black holes does not appear to be straightfor- 
ward. By applying Bohr's correspondence principle. Hod [28] has argued that the real part of 
the quasi-normal frequencies of gravitational waves ought to be given by the asymptotic expres- 
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sion (c/. eq. (1)) 

3?a; = Tffln3 + mQ (3) 

where m is the azimuthal eigenvalue of the wave and fl is the angular velocity of the horizon. On 
the other hand, numerical results have been obtained [10] which appear to contradict the above 
assertion, suggesting instead an asymptotic expression independent of the temperature, 

= mQ (4) 

Here, we present an analytical solution to the wave (Teukolsky [29]) equation which is valid 
for asymptotic modes bounded from above by 1/a, where a is the angular momentum per unit 
mass of the Kerr black hole. Thus, our calculation is vahd in the small-a limit (a <S 1), which 
includes the Schwarzschild case (a = 0). Our results confirm Hod's expression (3), but do not 
necessarily contradict the numerical result (4). The latter may well be valid in the asymptotic 
regime 1/a < uj. In the Schwarzschild limit {a — 0), the range of frequencies considered in our 
calculation extends to infinity and our expression, which generalizes (3), reduces to the expected 
form [5] generalizing (1) to arbitrary spin of the wave. 

The metric of a four-dimensional Kerr black hole may be written as 

,2 / 2Mr\ ,2 4:Mar sin^e , E ,0 v. ,.2 
ds^^ - 1 — df + dtd(f) + — dr^ + T.d9^ 

.2^/2 2 2MaVsin2^\ 
+ sm^ e\r^ + ^ dcf)^ (5) 



E 
where 

E = + cos^ 9 , A = - 2Mr + (6) 

M is the mass of the black hole and we have set Newton's constant G —1. The roots of A are 
given by 

r± = M ± VM2 - a2 (7) 
the larger being the radius of the horizon (r/j = r+). The black hole is rotating with frequency 



and its Hawking temperature is 



Q = (8) 
2Mr+ ^ ' 
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Small perturbations are governed by the Teukolsky wave equation [29-31] 



—2s — — — r — la cos 6 — — 

V A ) dt 

1 d ( . (9^\ fa{r-M) icos^\ <9^ . 2 2 



sinOan'^'wj-^^l^^+WflJ a^ + (^^cot^''-«)* = (10) 
where s = 0, —1, —2 for scalar, electromagnetic and gravitational wave, respectively. 
Separating variables 

* = e-^'"V"*'^5(^)/(r) (11) 
we obtain an angular and a radial equation, respectively, 

-^{smesy+( aWcos'^e- - 2awscose - '^H^f^ _ ^^cot^^^ S = -{A+s)S (12) 

sm6' \ sin^6' sin^6' / 

i^(A^+V')' + Vir)f = (A + aV)/ (13) 

where the potential is given by 

, (r'^ + a^fu"^ - 4:aMrujm + a^m^ + 2ia(r - M)ms -2iM(r'^ - 0^)0:3 ^ , 
V{r) = ^ ^ V 2irus (14) 

and A is an eigenvalue to be determined by solving the angular eq. (12). Before attempting to 
solve these equations, let us simplify them by placing the horizon at r = 1. This is accomplished 
by setting 

2M = 1 + (15) 
Then the roots of A (eq. (7)) are given by 

r- — , r+ — 1 (16) 

respectively. Let us try to solve the two wave equations (12) and (13) by expanding in a. We 
shall keep terms up to o(a). Moreover, we shall assume that the frequency u; is large but bounded 
from above by the large parameter 1/a, 

u < 1/a (17) 



4 



Thus, uj is not in the asymptotic regime but rather in an intermediate range. This range of 
frequencies includes the asymptotic regime in the Schwarzschild hmit a — > 0. Our calculation is 
valid in the limit of small a (a <^ 1). 

The solutions of the angular equation (12) to lowest order are spin- weighted spherical har- 
monics [29] and the separation constant (eigenvalue) A is 

^ = £(£ + 1) - s(s + 1) + o{auj) (18) 

It is convenient to express the radial eq. (13) in terms of a "tortoise coordinate." To define it, 
observe that near the horizon (r — > 1), the wavefunction behaves as 

/(r)~(r-l)^ , X = i{uj-am) + o{l/uj) (19) 

On the other hand, at infinity (r — > oo), we have 

/(r) ~ e'^^ (20) 

By introducing the variable ( "tortoise coordinate" ) 

z — ujr -\- {uj — am) ln(r — 1) (21) 

we may express the boundary conditions simply as 

f{z) ~ e^'' , z^±oo (22) 

Moreover, the monodromy for the function e'^^f{z) (which approaches a constant as r ^ oo (i.e., 
z —>■ +oo)) around the singular point r = 1 to this order (o(a)) is easily deduced from (19), 

M(l) = e^'^^'^-"'"^ (23) 

This will serve as the definition of the boundary condition at the horizon [5]. The contour 
surrounding the singularity r = 1 can be deformed in the complex r-plane so that it either 
lies beyond the horizon (3?r < 1) or at infinity (r — >■ oo). Then the monodromy only gets a 
contribution from the segment lying beyond the horizon. To express the radial equation in terms 
of the tortoise coordinate, define 



where Aq = r(r — 1) (note from eq. (6) that A = Aq + o(a^)). Inverting (21), 



r^\l-^ + o{l/u;) (25) 

we obtain the radial equation (13) to lowest order in Xj ^Juj in terms of i?, 

d^R (_ 3is 4-s^-Aiams] „ „ 

1 + — + -— } R^O (26) 



to be solved along the entire real axis. This is Whittaker's equation [32] . Two linearly indepen- 
dent solutions are (setting x — 2iz) 

M^,±i,{x) = e-^/^x^'^+^/^M^ ±1^- K,l±2i^,x) (27) 

where 

3s 2 g(g + ^iam) 
''^T ' ^ = ^6 (28) 

and M is Rummer's function (also called $). For our purposes, the following linear combination 

(Whittaker's function) will be useful: 

due to its clean asymptotic behavior, 

W^,^{x)r.e-^/' X- {l + o{l/x)) (30) 

as |x| oo. 

To compute the monodromy around the singularity r = 1, we shall deform the contour so that 

it gets mapped onto the real axis in the z-plane. Near the singularity 2; = 0, we have 2; ~ — ^ 

(eq. (25)). We shall choose our contour in the complex r plane so that near r — 0, the positive 

real axis and the negative real axis in the z-plane are mapped onto 

argo; 3% argo; 
argr = 7r — , argr = ^ — (31) 

in the r-plane. These two segments form a 7r/2 angle (independent of argo;). To avoid the r = 
singularity, we shall go around an arc of angle 37r/2. This translates to an angle Stt around z — 
in the 2;-plane [5]. To implement this analytic continuation, observe 

M^^±f,{e^^'x) = e3'^*(±'^+^/2) e-^/2^±/^+V2M(i ± + k,1 ±2i^,x) 

= -ie^'-'^M^,,±^(x) (32) 



where we used 

M{a,b,-x) = e^''M{b-a,b,x) (33) 

Therefore, 

For the asymptotic behavior, we need 

Jx) = ^ , e-^-'^W^^ Je'^x) + ^/^ ^ e-'"(i+'^+''%_, Jx) (35) 
As — > oo, we obtain 

p^i^^^^^e e ( ^) e e x (3bj 

and so 

W^^^^ie^'^'x) ~ Ae'/'^x^ + -Be-^/'a;-'^ (37) 

where 

« = - c - (3S) 

Using (28), after some algebra we obtain from (38) 

^= -(l + 2cos7rs) + o(a^) (40) 

where we also used the identities 

r(i - x)r(x) = , r(i + x)r(i - x) = (4i) 

Eq. (40) gives the correct Schwarzschild limit [5]. Notice that there are no o(a) corrections. 
Using (23), we obtain for the monodromy around r = 1, 

^gMc-ma) (42) 

therefore, 

3?a; = ln(l + 2 cos tts) + ma + o(a^) (43) 
47r 



which is in agreement with Hod's formula (3) in the case of gravitational waves (s = —2) and 
in the small-a limit (in which fl ^ a; also note Th ~ ^ in our units). However, it should be 
emphasized that these are not asymptotic values of quasi-normal modes being bounded from 
above by 1/a. 

To summarize, we have obtained an analytical expression for asymptotic values of quasi- 
normal frequencies of Kerr black holes by solving the wave (Teukolsky [29]) equation pertur- 
batively. The zeroth-order approximation is Whittaker's equation. We apphed the monodromy 
argument of Motl and Neitzke [5] to the zeroth-order solution and arrived at an explicit expression 
for arbitrary spin of the wave (eq. (43)). This result is in agreement with Hod's suggestion (3) 
based on Bohr's correspondence principle in the case of gravitational waves and in the small-a 
hmit. Eq. (43) is only apphcable to modes bounded from above by 1/a. It would be interesting 
to extend our results to general values of a and an unbounded asymptotic regime of quasi-normal 
modes, 1/a < lu. Such an extension is far from straightforward. 



Acknowledgments 

We wish to thank J. D. Bekenstein, E. Berti, V. Cardoso, K. Kokkotas, L. Motl and A. Neitzke 
for useful discussions. G. S. is supported by the US Department of Energy under grant DE- 
FG05-91ER40627. 



8 



References 

[1] O. Dreyer, Phys. Rev. Lett. 90 (2003) 081301; gc-qr/0211076. 

[2] G. Kunstatter, Phys. Rev. Lett. 90 (2003) 161301; gr-qc/0212014. 

[3] L. Motl, Adv. Theor. Math. Phys. 6 (2003) 1135; gr-qc/0212096. 

[4] A. Corichi, Phys. Rev. D67 (2003) 087502; gr-qc/0212126. 

[5] L. Motl and A. Neitzke, Adv. Theor. Math. Phys. 7 (2003) 2; hep-th/0301173. 

[6] E. Berti and K. D. Kokkotas, hep-th/0303029. 

[7] R. A. Konoplya, Phys. Rev. D68 (2003) 024018; gr-qc/0303052. 

[8] A. Maassen van den Brink, gr-qc/0303095. 

[9] A. P. Polychronakos, hep-th/0304135. 

[10] E. Berti, V. Cardoso, K. Kokkotas and H. Onozawa, hep-th/0307013. 

[11] V. Cardoso, O. J. C. Bias and J. P. S. Lemos, Phys. Rev. D67 (2003) 064026; 
hep-th/0212168. 

[12] S. Chandrasekhar and S. Detweiler, Proc. R. Soc. London, Ser. A 344 (1975) 441. 

[13] E. W. Leaver, Proc. R. Soc. London, Ser. A 402 (1985) 285. 

[14] H. P. NoUert, Phys. Rev. D47 (1993) 5253. 

[15] N. Andersson, Class. Quant. Grav. 10 (1993) L61. 

[16] A. Bachelot and A. Motet-Bachelot, Annales Poincare Phys. Theor. 59 (1993) 3. 

[17] S. Hod, Phys. Rev. Lett. 81 (1998) 4293; gr-qc/9812002. 

[18] J. F. Barbero G., Phys. Rev. D51 (1995) 5507; gr-qc/9410014. 

[19] G. Immirzi, Nucl. Phys. Proc. Suppl. 57 (1997) 65; gr-qc/9701052. 

9 



[20] C. Rovelli and L. Smolin, Nucl. Phys. B442 (1995) 593; gr-qc/9411005. 

[21] A. Ashtekar and J. Lewandowski, Class. Quant. Grav. 14 (1997) A55; gr-qc/9602046. 

[22] C. Rovelli, Living Rev. Rel. 1 (1998) 1; gr-qc/9710008. 

[23] T. Thiemann, gr-qc/0 110034. 

[24] C. Rovelli and R Upadhya, gr-qc/9806079. 

[25] J. D. Bekenstein, Lett. Nuovo Cimento 11 (1974) 467. 

[26] J. D. Bekenstein and V. F. Mukhanov, Phys. Lett. B360 (1995) 7. 

[27] S. Musiri and G. Siopsis, hep-th/0308168. 

[28] S. Hod, Phys. Rev. D67 (2003) 081501; gr-qc/0301122. 

[29] S. A. Teukolsky, Phys. Rev. Lett. 29 (1972) 1114. 

[30] S. Mano, H. Suzuki and E. Takasugi, Prog. Theor. Phys. 95 (1996) 1079; gr-qc/9603020. 

[31] M. Sasaki and H. Tagoshi, gr-qc/0306120. 

[32] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series and Products, Academic Press 
(1994). 



10 



